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THE EQUILIBRIUM THERMODYNAMIC PROPERTIES 
OF A HIGH-TEMPERATURE NITROGEN— ALKALI- METAL 

VAPOR MIXTURE 

By Stephen K. Park 
Langley Research Center 

SUMMARY 

A study has been made of the equilibrium thermodynamic properties of nitrogen, 
seeded with alkali- metal vapor, at high temperatures (1000° K to 20 000° K), moderate 
pressures (0.003 to 3 atmospheres), and 0- to 5-percent seeding. It was found that 
because of the low ionization potentials of the seed materials, their ionization could repre- 
sent a significant energy sink in regions where nitrogen is essentially chemically inactive. 
In addition, the increased number density of electrons resulting from ionization of the 
seed material causes a delay, until a higher temperature, in the ionization of nitrogen. 
These two effects tend to produce, for a given seed material, rather well-defined 
temperature-pressure regions in which even small degrees of seeding produce fairly 
large variations in some thermodynamic properties. In other regions, if proper allow- 
ances are made for the change in the average atomic weight of the mixture (and conse- 
quently, a change in the gas constant of the mixture) caused by the mass fraction of seed 
material, the effect of seeding is found to be very small even for as much as 5-percent 
seeding. 

In general, it was found that the effect of seeding on the thermodynamic properties 
of the mixture tends to become more pronounced at lower pressures. Also, as might be 
expected, the effect usually becomes more pronounced as the seeding percentage is 
increased. In addition, certain properties (the specific heats and, to a lesser degree, 
enthalpy and internal energy) are significantly more sensitive to seeding than others 
(entropy, speed of sound, departure coefficient, for example). Although the most pro- 
nounced effects are associated with cesium seeding, effects caused by potassium and 
sodium seeding are very similar and only slightly less pronounced. 

A computational procedure is presented whereby the thermodynamic properties of 
this mixture may be computed for arbitrary degrees of seeding. This procedure is 
readily adaptable to high-speed computer use and should be useful as a subroutine in 
connection with step-by-step integration solutions of high-enthalpy gas flow problems. 
Included in this procedure is an approximate general method of treating electronic exci- 
tation and variable principal quantum number cutoff. This method permits the calculation 



of excitation to high principal quantum numbers without a large increase in computation 
time. 


INTRODUCTION 


In recent years there has been a considerable amount of research with high- 
temperature plasma devices, such as plasma accelerators and magnetohydrodynamics 
(MHD) generators. Some of these utilize a nitrogen — alkali-metal vapor mixture as the 
working fluid because the addition of a small amount of alkali- metal vapor can result in 
a large increase in the electrical conductivity of the gas. Usually, an equilibrium solution 
of the equations governing the operation and performance of these devices requires a 
high-speed computer program using the thermodynamic properties of the gas mixture in 
connection with a step-by-step integration procedure. In addition, some of these devices 
(see, for example, ref. 1) produce gas enthalpies corresponding to those associated with 
reentry velocities and can be used in this area to study heat transfer and related problems. 
Although there are available a number of calculations of the thermodynamic properties of 
high'-temperature nitrogen (refs. 2 to 7, for example), little data on the thermodynamic 
properties of a high-temperature nitrogen — alkali-metal vapor mixture are available. It 
is therefore desirable to develop expressions that will accurately describe the properties 
of this mixture and that will be in a form easily adapted for use with a computer program. 

In this paper expressions for the thermodynamic properties of a high-temperature 
(up to 20 000° K) nitrogen — alkali-metal vapor mixture are developed by using a partition 
function analysis. This method is similar to, for example, that employed by Drellishak, 
et al. (ref. 2) for nitrogen, oxygen, and argon and by Ahtye and Peng for nitrogen (ref. 3). 
The reactions considered are the dissociation and single ionization of nitrogen and the 
single ionization of the alkali- metal vapor. 


Expressions are obtained for the following conventional properties: entropy, 
enthalpy, internal energy, specific heats at constant pressure and volume, departure 
coefficient, density, and speed of sound. In addition, expressions are obtained for the 
logarithmic partial derivatives of the departure coefficient with respect to temperature 


and pressure R 


[i + lM'il 

and R 

1 . P " 
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, quantities which are fundamental to 


real- gas flow problems and appear as a consequence of the dissociation and ionization 


reactions. 


A treatment of the electronic excitation partition function and its first and second 
derivatives is presented. Included in this treatment are an estimation and compilation of 
unobserved energy levels and the development of a simple closed-form expression to 
estimate the high principal quantum number contribution to the partition function. A 
method, which is a simplification of the method suggested in reference 8, of treating the 
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problem of a variable principal quantum number cutoff is discussed. This method permits 
the easy handling of a variable cutoff without the usually associated large increase in 
computation time. 

An outline is presented herein of a computer program which calculates the thermo- 
dynamic properties of a nitrogen — alkali-metal mixture as a function of temperature, 
pressure, and degree of seeding. This program, as presented, may be conveniently 
adapted for use as a subroutine in connection with computer programs used to solve real- 
gas high- enthalpy flow problems. 

Because they are commonly used as seed materials a rather extensive study has 
been made of the effect of cesium, potassium, and sodium seeding on the thermodynamic 
properties. It is shown that in certain temperature-pressure regions, even the addition 
of small amounts of seed material (1 to 2 percent) may cause large variations (20 to 
40 percent) in some of these properties. Results are presented for a temperature 
range of 1000° K to 20 000° K, at pressures of 0.003, 0.01, 0.03, 0.1, 0.3, 1, and 
3 atmospheres, and for degrees of seeding from 0 to 5 percent. 

SYMBOLS 


a speed of sound, m/s 

a number of parent configurations 

Cp specific heat at constant pressure, J/kg-°K 

c v specific heat at constant volume, J/kg-°K 

e electron 

E(i) internal energy of ith species, J/m3 

E total internal energy, J/kg 

f(T,n) any general function of temperature and principal quantum number cutoff 

statistical weight of vth state 
h Planck's constant, J-sec 

H total enthalpy, J/kg 
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Boltzmann constant, J/°K 


K w = m x /m N 




degree of seeding (see eq. (7)) 


K 1= Kc 


1 + K r 


K r (Y) 

m i 

m 

n i 

n 

n t 

n o 

n 2 

N 

N+ 

N(n) 

Pi 

P 


dimensionless equilibrium constant 
mass of a particle of species i, kg 
principal quantum number 
number density of species i, particle/m 3 
principal quantum number cutoff 
total number density, particles/m 3 
reference quantum number 
nitrogen molecule 

nitrogen atom, also nitrogen generally 
singly ionized nitrogen atom 

last term in excitation partition function summation 
partial pressure of species i, N/m 2 
total pressure, N/m2 


p Q atmospheric pressure, N/m2 

P = P/P 0 > atm 

q effective charge 
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Qj (i) partition function of ith species 

Q(i) total partition function of ith species 

R mixture gas constant, J/kg-°K 

S(i) entropy of species i, J/m3-°K 

S total entropy, J/kg-°K 

T temperature, °K 

V volume, kT/p, m3 

X alkali-metal atom, also an alkali metal generally 

X + singly ionized alkali- metal atom 

Y type of gas 

Z departure coefficient, p/pRT 

Z T = R f + |(^) p ]’ J ^-°K 

z P = R f-#(f)^ J/k e-° K 

a r degree of reaction (r = 1,2,3) 

Po = f + 6 rot( N 2) + e v( N 2) + e e( N 2) + 2K C | + e e (X)J 

Pr = T art l0 ^e K r(Y)] 

y 

A 
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L 


e j(i)= T (£l°g e Q)< 1 )) v 


0jj(Y) characteristic temperature of dissociation, °K 

q ^ _ Reference energy of ith species o K 
° k 

0 r (Y) characteristic temperature of reaction r, °K 

0„ characteristic temperature of energy level v, °K 

0 V characteristic temperature of vibration, °K 

0 rot characteristic temperature of rotation, °K 

0jj Rydberg constant, °K 



v summation index of excited states 

p total density, kg/m^ 



o> general seeded thermodynamic variable 

c o 0 general unseeded thermodynamic variable 

Subscripts : 


i species (N 2 , N, N + , X, X + , e) 

j energy forms: translational t, rotational rot, vibrational v, and electronic 

excitation e 

r reaction number (=1,2,3) 
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EQUATION OF STATE 


This paper considers a mixture of nitrogen and an alkali vapor in a pressure- 
temperature region (see fig. 1) where the principal reactions r are 


N 2 r N + N 

(r = 1) 

(1) 

x:x + + e 

(r = 2) 

(2) 

N H N 4- + e 

(r = 3) 

(3) 


and X denotes any alkali- metal atom. It is assumed that each species i (i being 
N 2 , N, N + , X, X + , e) obeys an ideal gas equation 


Pi = njkT 

where Pi is the partial pressure and ^ the number density of species i. (The 
validity of this assumption is discussed in ref. 8 where it is shown that for the range of 
temperature and pressure considered herein, the assumption is well justified.) The 
total gas pressure p then becomes 


P = S Pi = S ihkT (4) 

i i 

The total gas mass density p is the sum of the species mass densities so that 

p = S n^i 
i 

where mi is the mass of a particle of species i. If it is noted that 
1 m N2 = m N ~ m N +, m x « m x +, n e m e « p, and K w is defined K w = jJjX, 
equation (4) may be written as 


p = pRTZ (5) 

where 

r = h (6a) 

2nijj(l + K w Kc) 
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(6b) 


Z = 


2 ? n i 

2"N 2 + "N + ”N+ 


The quantity K c is the degree of seeding and may be defined in terms of number densi- 
ties as 


K c - 


n X + n X + 


2n N 2 + n N + n N+ 


(0 g K c < oo) (7) 


Some authors may prefer to define a slightly different degree of seeding K c as 


K c ls 


n X + n X+ 


(O g Kc 1 g l) 


( 2n N 2 + n N + n N + ) + ( n X + n X + ) 
but for small degrees of seeding, the choice of definition is unimportant since 


K c l = — ~ K c 

C K 0 + 1 c 


In general, the alkali-metal vapors (lithium, sodium, potassium, rubidium, and 
cesium) have a relatively high vapor temperature and a low characteristic temperature 
of dissociation. As a result, even at temperatures close to the liquid- vapor saturation 
line, the vapor is essentially completely dissociated. (See, for example, cesium, ref. 9.) 
In addition, the characteristic temperature of second ionization is on the order of, or 
much greater than, the corresponding characteristic temperature of second ioniza- 
tion for nitrogen. Therefore in the pressure-temperature region shown in figure 1, the 
predominant reaction involving the alkali metals is X ~ XT 1- + e. Values of K w and the 
characteristic temperatures of dissociation and first and second ionization are presented 
in table I. 


PARTITION FUNCTIONS 

The theory of partition functions as a consequence of Maxwell- Boltzmann statistics 
and its application to equilibrium computations for gas mixtures have been discussed in 
many references (for example, refs. 8 and 10 to 14). For this reason, partition functions 
and the related number density equations are introduced with a minimum of theoretical 
explanation. 
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The partition functions Qj(i) for species i possessing energy of the form j are 


Qt(i) = 




N(n) 

Qe(i) - ^ Sv e T 

v=l 


(8a) 


(8b) 


QrotW = ( 8c ) 

^rot 

( -M ’ 1 

Q v (N2) = V 1 ' e T / (8d) 

where Q ro t( N 2) is the classical limit of the rotation partition function and Q v (N 2) is 
the harmonic oscillator partition function. The energy forms are translation, electronic 
excitation, rotation, and vibration and the 6 values are characteristic temperatures of 
electronic excitation ^ Energy ^ level — rotation, and vibration. (See table I.) The 
only species containing energy in rotation and vibration is, of course, N2. The total 
partition function Q(i) for species i is 

Q(i) = TT Qi(i) 

j 

Partition functions Q ro t( N 2) and Q V ( N 2) are written in a rather simple form 
which ignores such real gas effects as anharmonic vibration and vibration-rotation cou- 
pling. In reference 15 it is shown that for equilibrium computations of gas thermody- 
namic properties in the pressure-temperature region considered herein, these effects are 
small (generally 1 percent or less) and they were therefore ignored. 

Since in the partition function Q e (i) the quantity N(n) is a function of the cutoff 
principal quantum number n, which is, for a given species i, dependent generally on 
temperature and electron and ion densities (ref. 16), Q e (i) is not a function of tempera- 
ture alone. However, a discussion of the excitation partition function is presented in 
appendix A where it is shown that 

(Jr >°Se Qe(») v * (w >°Se <5e<‘>)„ 
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and that 


_9_ 

9T 


T2(^ log e Qe(i))l - log e QeG))* 


n 


The evaluation of the thermodynamic properties of the gas mixture involves evalua- 
tion of the terms 

n 


T^loge Q(i)j v = T £ Z log e Qj(i) 


V 


and 


9_ 
1 9T 


T 2 ^ log e Qft) 


-±)t> 


Id 

\9T 


T2 fe? log e Q j (i) ' 


V 


The differentiation may be performed to obtain 


■pjr log e Q(i) - g + e eW + e v® + e rot^ 


( 9 ) 


where 


e e (i) 



T 


(i * e) 


e v (i)=^V eTV -V (i = N 2) 

e v (i) = 0 (i * N 2 ) 

«**« = 1 {i = N 2) 

e rot (i) = 0 (i* N 2 ) 


and to obtain 




= | + cr rot (i) + CT v (i) + CT e (i) 


( 10 ) 
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where 


°rotW - e rotW 


o- e (i) = 



v=l 



Q e (i) 



°v(i) = e v 2 (N 2 )e T 


(i * e) 


(i = N 2 ) 


a v (i) = 0 


(i * N 2 ) 


As is shown in appendix A, the quantities Q e (i), ee(i), and a e (i) are functions of 
both temperature and principal quantum number cutoff n. Shown in figures 2 to 4, 5 to 7, 
and 8 to 10 are Qe(i), e e (i), and <r e (i), respectively, plotted against temperature for 
various values of n. For the alkali metals, as is suggested by equations (A3) to (A6), 
for higher values of n, these quantities become functions of t/ 02( x ) only anc * therefore 
for each quantity only one plot is needed to describe cesium, potassium, and sodium (and 
presumably also lithium and rubidium). 


EQUILIBRIUM CONSTANTS 


For each of the equilibrium reactions 1, 2, and 3, a dimensionless equilibrium 
constant may be defined as times a function of the number densities of the species 

involved. This function may be related to the ratio of the respective total partition func- 
tions to obtain 


Pq. 

kT 


K X (N) 


o M N) 

n N 2 _ Q 2 ( n ) e “ T 
n N2 VQ(N2) 


(11a) 


£o_ 

kT 


K 2 (X) 


e 2 (x) 

n X-^e = Q(X + ) Q(e) e “ ~T~ 
n X VQ(X) 


Pq. 

kT 


Kg(N) 


0 3 ( n ) 

"N^e Q(N + ) Q(e) " T 
n N VQ(N) 


(lib) 


(11c) 
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where K r (Y) is the dimensionless equilibrium constant, and 0 r (Y) is the characteris- 
tic temperature, for reaction r involving species of gas Y = (nitrogen, N or al k ali 
metal, X). The quantity 0 r (Y) is defined in terms of the net energy change Ae r (Y) 


. v Ae r( y ) 

for reaction r, as 0 r (Y) s — - . 

may be expressed as 


/irm N k\ 3 / 2 

Kl(N)= h^ 


By using equations (8a) to (8d), the K r (Y) values 

0l(N) 


/k_\/^rot\(i _ e ~ t)t 5 / 2 
v °/\ T / ' ' Q e (N 2 ) 


K 2 (X) = 2(t^f /2 (±y/2 Sem e - T 
\ h 2 / [PoJ Q e ( x ) 


e 2 (x) 


K 3(N ,. 2 (^) 3/2 (i) T V2^) e - 9 f 

where Q e (e) = 2 (e e (e) = cr e (e) = o). Shown in figure 11 are values of K r (Y) plotted 
against temperature for various pressures and K c = 0.01. 

Some important derivatives which will be needed later are given by 

Pi = T ^r[log e Kj(N)] v = | - e rot (N 2 ) - e v (N 2 ) - e e (N 2 ) + 2e e (N) + (12a) 


H s T If |log e K 2 (X)] v = | - e e (X) + e e (X+) + 9 -l~ (12b) 

P 3 = T ^|loge K 3 (N)] = | - e e (N) + e e (N+) + ^5. (12c) 

and 

H s w[( T ^l)]p = | ' a rot(N 2 ) - a v (N 2 ) - a e (N 2 ) + 2a e (N) (13a) 

V-2 £ w[( T ' 32 )]p = I " a e( X ) + °e0?) (13b) 

p 3 h ^[(T/3 3 )] p = | - °e(N) + Oe(N+) (13c) 
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DETERMINATION OF SPECIES CONCENTRATIONS 


In general, thermodynamic properties can be expressed as functions of any two of 
the three variables - temperature, pressure, and density. The equilibrium properties of 
a particular species are conventionally described by the temperature and it is primarily 
through the solution of the number density equations (Ha) to (11c) that the pressure or 
density dependence is introduced. Since it is the more conventional property and since 
in most applications it is more easily measured, pressure rather than density was chosen 
herein as the second independent variable. 

Equations (11a) to (11c) and (7) together with expressions for charge neutrality 

n e = n X"*" "** n N + 

and species conservation 

form a set of six equations in the six variables n^. If for the reactions r = 1,2,3, 
respectively, a degree of reaction a? r (o < a r < l) is defined as 

n N 

a, = „ 

1 2n N 2 + n N 

01 2 n x + n x + 


s n N + 

"3 n N + n N + 

the set of six equations can be reduced to three equations 


Kj(N) 2a x 

(1) 

«1 - 0-^3 

P 1 - o'! 

w 

1 - « 3 (l - «i)_ 


K 2< X > «2 A\L , a l 0! 3 

P 1 - a 2 \ ZJ c 2 l-a 3 (l- ai ) 


(14a) 


(14b) 
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(14c) 


k 3 (n) 

P 


Oi< 


1 - agVZ, 


K c“2 + 


a l a 3 


1 - a 3 ( 1 - o-j) 


p 

where P = — and 

Po 


It is shown in appendix B and illustrated in figure 1 that reactions 1 and 2 may occur 
in overlapping temperature-pressure regions, but that reaction 3 becomes significant 
(oi % ± Oj only in temperature-pressure regions where reactions 1 and 2 are essentially 
complete (ai ~ a 2 ~ l). Therefore aj, a 2 , and 0-3 may be expressed to a high degree 
of accuracy as 


01 + 2 q' 1 q'« 


Z = 1 + 2K C (1 + a 2 ) + 1 7 1 6 

1 - “ si 1 - “ l ) 


“1 = 


1 + 4 


1 - 1/2 


(1 + A)K X (N) 


(15) 


Q!o = 


f 1 + “l ] k 2 (x) 
V 4 ) KcP 


(1+A) 


1 + 


8K C P 


(l + ai )(i + A)K 2 (X) 


- 1 


(16) 


K C (P + 2K 3 (N)) 

3 " °(p + k 3 (n)) K 


a o = 


1 + 


4(1 + 2K C )(P + K 3 (N))K 3 (N) 
K C 2 [P + 2K 3 (N)] 2 


(17) 


1 + a 


Note that in this form, equations (15) and (16) are still coupled, 


where A = 2K„ 

c 1 + a 

but, since K c < A < 4K C and usually K c « 1, equations (15) and (16) are very insensi- 
tive to A and an iteration solution will therefore rapidly converge. Since the K r (Y) 
terms are functions of pressure (through the principal quantum number cutoff, n), an 
iteration solution on n is also necessary. However, since the solution of equations (15), 
(16), and (17) is very weakly coupled to n (through n e ” 1/ /, this iteration procedure 
also rapidly converges. It follows that for a given temperature, pressure, and degree of 
seeding, o-j, a 2 , and a 3 are uniquely determined and the quantities — 


ni 


may be obtained from the following relationships : 


N 0 + n N + n N + 
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n i 

2n N 2 + n N + n N+ 


2 

(i = N 2 )^ 

“1 - “l“3 

(i = N) 

00 

Q 

t-H 

Q 

(i = N+) 

K c (l - a 2 ) 

S' 

ll 

K-c a 2 

II 

3 

c a 2 + 01 ^ci 3 

(i = e)J 


The number density ni of species i is then given by 


and 


n- = n i 2 P 

2n N 2 + n N + n N+ Z kT 

Z = 1 + Q!j + 2K c (l + 0i2 ) + 2Q!jO!g 


( 18 ) 


(19) 


It is interesting to note the effect of electronic excitation on the solution for o; r . 
Equations (15) to (17) were solved (K c = 0.01) by using K r (Y) values which were insensi- 
tive to electronic excitation ^all Q e (i) values were set equal to their ground- state degen- 
eracy q^o) anc * va l ues a r obtained are shown in figure 12 for pressures of 0.01 
and 1 atmospheres. Note that for small values of a r , the effect is negligible but that as 
o! r approaches 1, the effect becomes apparent. This effect is caused by excitation 
becoming an appreciable energy sink at high temperatures. 

The first derivatives of a r are obtained from equations (B6) to (B8) in appendix B. 
To obtain the first partial derivatives of <23, equation (B8) may be differentiated directly 
to obtain 



(20a) 


T_ 



p 


03 


(20b) 


15 



but, because of the coupling, the first derivatives of and cn 2 must be obtained from 

the differentiation of equations (B6) and (B7) and the simultaneous solution of the equations 



/da 


a 


9p 


' B 0t 

V“ap" 




Cl 


'T,p 


^b 

a P 


where r = a = 1 when r = b = 2 and r = a = 2 when r = b = 1. By performing the 
indicated operations and simplifying, it is possible to obtain 


(tt)p = ^t ^1 ' 2 ^ 


(20c) 


\9p. 


= 'P^1 


(20d) 


(dot o\ j3o , 

[it) d = T ^2 + t( 0 i ' ^2) 


(20e) 


tel = ' P ^2 


(20f) 


^The equations for to cj)^ are given in appendix B.) 

Note that the indicated differentiation leading to equations (20a) to (20f) has been 
performed under the assumption that for 0 < a r < 1, the quantities Q e (i) (and therefore 

K r (Y)) are weak functions of pressure ~~~j « lj- This assumption has been 

verified numerically as would be expected since generally it is only for a r ~ 1 that 
K r (Y) becomes pressure dependent. (See figs. 11 and 12.) 

INTERNAL ENERGY, ENTHALPY, AND ENTROPY 


In terms of the total partition function, the internal energy per unit volume for 
species i is given by 
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where 


E(i) = nikT 2 f^ log Q(i)l + ni k0 o (i) 
y J- J V 


0 o (i) = 0 (i = N 2 ,tf) 

= \ 6 1< N > (i = N) 

0 o (i) = 0 2 (X) (i = X+) 

9 0 (i) = \ 0 l( N ) + «s(N) tt = N 4 ) 

0 o ( i) = 0 (i = e) 

(values of 0 r (Y) are given in table I.) 


Note that the reference energy is taken to be zero for the nitrogen molecules N 2 
and alkali- metal atoms X. It was felt that in most engineering applications, the seed 
material would be introduced into the system at high enough temperatures to assure the 
material would be in a monatomic vapor form. Therefore, the reference state for the 
enthalpy H and internal energy E was taken to be a gas mixture of N 2 and X. If 
it is desired to reference an enthalpy H' and internal energy E' to a gas mixture of 
N 2 and X 2 , it may be done with the following relation: 


H’ - H = K c R0 d (X) = E' - E 


The total gas energy per unit mass E is therefore 

S E(i) E E(i) 

jTj „ _i „ _i 

p E n^mj 

i 


which becomes 



where R is defined by equation (6a). 

The gas enthalpy H is simply determined by 


H = E+ — = E + RTZ 
P 


( 21 ) 


(22a) 


17 



With the use of equations (18), (11a) to (11c), and (12a) to (12c), H may be expressed as 


H = RT(|3 0 + a 1 ^ 1 + 2o' 1 a'3^3 + 2K c a 2$2) 


where 


ft o*| + e r ot< N 2> + + 6 e< N 2> + 2K c[| + «e(»] 


The entropy per unit volume S(i) of species i is 

S(i) = ink^loge ^ + 1^ + nikjT log e Q(i)J 
so that the total entropy per unit mass S is 


V 


S - 2R 1 2 n N rn N -- n ^( l0g e * [ T W l0 «e «W>] v} ♦ RZ 

i 2 


or from equations (21) and (22a) 

H 


S T + 2R X 2n N + n N + n N + 
i * 


log 


Q(i) 0 o (i) 


e Vm 


SPECIFIC HEATS AND SPEED OF SOUND 


The specific heat at constant pressure c p is defined as so that from 

equations (13a) to (13c) and (22), ^ 

c p ( da -\\ 

= Ho + a + 2q' 1 q' 3P3 + 2K c a 2 P2 + T ^l^afy 


+ 2T/3 3 


/9a, 

) 

\ /P 


a A-d^l + a 3[W L > 


+ 2Kcft 2 x(^) p 


where 


^o 


■[M 


Pp = k + CT rot^2) + a v(N2) + a e(N2) + 2K C + a e (X) 


(22b) 
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The specific heat at constant volume c v may be obtained from the relationship 
(ref. 13) 


T /9V\ 2 / 3V\~ _ tM 2 M 
’p - c v - - pV^ 9 T J p ^ g p L " p 2 \3T/ p \3p/ r 


2/a n \“l 
T 


( 23 ) 


so that 


c D - c v = z 


where 


Zip = R 


1 + 



= r T/W 

R V\9T 


Zp = R 


i - pfczN 

Z W/ T 


= - R 


P_/3V' 

V\9p 


and where the partial derivatives in equation (23) are obtained from equation (5). 
The speed of sound is given by 


i.2 = it M = yRTZ 
P\B log e P/g 


(24) 


where 


y - 


8 lQ g e p \ = V_r\ 

9 lo §e P/g c v\ZpJ 


(25) 


Equation (25) is easily obtained from the differentiation of equation (5), the definition 
of an isentropic process, and the identity 

ifm) 


Z \3T 


9P/ n 


as 


The quantities Z T and Z p can be obtained from equations (19) and (20a) to (20f) 

IT = 1 + zll 1 + 2a 3)^1^1^4 + 2of 1^3^3 + 2 k c^2^2^4] 


— 1 + ^ j(^ + + ^• a x < P^ + ^K c ^>2^* 4] 
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As previously mentioned, the equilibrium flow solution of the equations governing 
the performance of plasma accelerators and MHD generators is usually accomplished 
with a computer by using a step-by-step integration procedure. The choice of what ther- 
modynamic properties are required at each step in the integration is, of course, dependent 
upon the form in which the flow equations are written. If, for one-dimensional MHD 
accelerator (or generator) flow, these equations are written as suggested in reference 17, 
practically all the real gas effects are explicitly contained in the property 0 which is, in 
terms of the variables used herein, 

Zrrf 

c p 

This quantity (3 is more than a convenient arrangement of terms and is, in fact, the real 

Cp 

gas analog of the ideal gas term - 1- 


RESULTS AND DISCUSSION 


The addition of a small amount of alkali- metal seed material into a nitrogen gas 
produces several rather distinct effects. Because of the low ionization potential and 
small partial pressure of the seed material, a greatly increased number density of elec- 
trons can be realized in a low-enthalpy temperature-pressure region. Shown in figure 13 
is the resulting number density of electrons obtained by seeding with various amounts of 
cesium, and in figure 14 the cesium electron densities are compared with those obtained 
by sodium and potassium seeding. Much of the ionization of the alkali metals occurs in 
temperature-pressure regions where nitrogen is chemically inactive and, therefore, in 
these regions represents a significant energy sink. In addition, the increased number 
density of electrons resulting from this ionization causes a delay, until a higher tempera- 
ture, in the ionization of the nitrogen atoms. Also, the addition of a seed material with a 
mass fraction will cause a reduction in the mixture gas constant R (see eq. (6a)) 

and all quantities explicitly containing R. This effect, although important (especially for 
cesium with K w ~ 10), is easily taken into account and for this reason it is more advanta- 
geous to study the seeding-induced variation in some thermodynamic properties oj by 
actually considering the property o>/R. If proper allowances are made for changes in 
the gas constant, it is found that the specific heats (and to a lesser degree, enthalpy, inter- 
nal energy, and the parameter 0 ) are significantly more sensitive to seeding than others 
such as entropy, speed of sound, Z-j/R, y, and Z. (Zp/R is almost completely insen- 
sitive.) The seeding-induced variations in Cp/R, H/R, 0, Z^R, S/R, Z, and y 
are shown in figures 15 to 21 for cesium seeding. In each case the quantity shown plotted 
is co/o> 0 or (w/R)/(o)/R) 0 = (1 + K w K c )(w/o) 0 ) where a> 0 is the corresponding prop- 
erty for unseeded nitrogen. Generally, as might be expected, the variations induced by 
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sodium and potassium seeding are very similar but are less pronounced and occur at 
slightly higher temperatures. This effect is illustrated in figure 22 for the specific heat 
Cp/R at pressures of 0.003 and 0.3 atmosphere. Shown in figures 23 to 29 are plots of 
<u 0 = c p> H, /3, Z-p, S, Z, and y. From figures 15 to 29 and values of K w (see 
table I), values of co may be obtained for various degrees of seeding. 

As a further note of explanation it should be emphasized that in figures 15, 16, 18, 

and 19, the quantity plotted is (1 + K W K C )-^- and that to obtain the actual value of co/w Q , 

Ct>o 

the factor (1 + K W K C ) must be taken into account. For example (from fig. 15), at a pres- 
sure of 0.03 atmosphere and a value of K c = 0.05, the quantity |^p/( c p)^] (■*■ + Kw^c) is 
approximately 1.8 for T = 3500° K while at T = 7500° K this same quantity is only 
0.8. Since for cesium with K c = 0.05 the quantity (1 + K w Kc) * 1.47, the actual 
values of cp^cp^ are 1.22 and 0.54. 

As mentioned previously, the partition function analysis employed herein is similar 
to that used in references 2 and 3. There are, however, some significant differences. In 
reference 3, Ahtye and Peng assume that reactions r = 1 and r = 3 are independent 
and they use an approximate method of fairing to compute from one region of reaction to 
another. In addition, they assume a constant principal quantum number partition function 
cutoff of n = 5. In contrast, Drellishak, et al. (ref. 2) include: a variable cutoff (only up 
to n = 32, however), no assumptions of independence of reaction, lowering of the ioniza- 
tion potential, and a more exact model of the N2 molecules. It may be seen that in 
comprehensiveness, for unseeded nitrogen, these two analyses bracket the analysis used 
herein. The data in table n are presented so that a comparison may be made of these 
analyses. 

As mentioned in the section "Determination of Species Concentrations" the solution 
for a r requires an iteration procedure involving Q e (i). Because of the form of Q e (i) 
this iteration procedure would normally be very time consuming. However, in connection 
with the discussion in appendix A, this iteration may be rapidly completed, for a fixed 
temperature and pressure, by computing Q e (i)j and then simply correcting the value of 
Q e (i) at each iteration by computing Q e (i )2 (eq. (A4)). 

It has been mentioned previously (for example, ref. 15) that a constant cutoff of 
n = 5 is usually considered to be sufficiently accurate for the computation of thermody- 
namic properties. As a check on this assumption and by using a computer program simi- 
lar to that presented in appendix C, certain thermodynamic properties were computed by 
using a constant value of n = 5 (values for the species N were taken from ref. 15) and 
these properties were compared with those obtained by the more accurate variable cutoff 
method (appendix A). These comparisons are presented in table m for unseeded nitrogen 
and also for 5-percent cesium seeding (n = 5 was also used for the cesium atoms although 
a larger n would have been more appropriate). As may be seen, fairly large differences 
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may occur. Inclusion of a reduced ionization potential would slightly reduce this effect 
but because of the present uncertainty in what the actual ionization potential reduction 
should be, it was not included in the calculations. 

CONCLUDING REMARKS 

Equations describing the thermodynamic properties of a high-enthalpy nitrogen, 
alkali- metal vapor mixture may be expressed in relatively simple form. This is true 
even if such non- ideal gas effects as electronic excitation and the associated variable 
quantum number partition function cutoff are included. The resulting expressions lend 
themselves well for use as a computer program subroutine (an outline of which is pre- 
sented herein) for supplying thermodynamic data in connection with step-by-step integra- 
tion procedures. 

In certain high-temperature devices using nitrogen seeded with an alkali-metal 
vapor, it is usually assumed that because of a lack of data on nitrogen — alkali- metal mix- 
tures and since the degree of seeding is small, the thermodynamic properties of this mix- 
ture are essentially those of nitrogen. It is found, however, that in certain temperature- 
pressure regions, 1/2- to 1-percent cesium seeding, for example, can change certain 
properties (such as the specific heats) by as much as 10 to 20 percent. In general, these 
effects become more pronounced at low pressures and/or high degrees of seeding. 
Although the largest effects are associated with cesium seeding, these effects become 
only slightly smaller for potassium and sodium seeding. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 28, 1967, 

129-01-05-04-23. 
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APPENDIX A 




ELECTRONIC EXCITATION PARTITION FUNCTIONS 


There are several problems encountered in the evaluation of the electronic excita- 
tion partition function Qe(i) and its derivatives e e (i) and cr e (i). The first of these is 
that in references 18 to 20, particularly at high principal quantum numbers, many pre- 
dicted energy levels have not been observed and are therefore not tabulated. For the 
alkali-metal atoms and also the nitrogen atom, failure to include these unobserved energy 
levels can generate, at high temperature, significant differences in Qg(i) and its deriv- 
atives. For the ions X + and N + , however, this effect is not important at temperatures 
below at least 20 000° K. To estimate the missing energy levels for the alkali-metal 
atoms it was assumed that for high principal quantum numbers and even at relatively low 
principal quantum numbers but high azimuthal quantum numbers, an energy state 6 V 
with principal quantum number m was hydrogen-like and could be, with reasonable 
accuracy, estimated by 

e p = e 2 (x) - (A1) 


where = 157 890° K. For the nitrogen atom, missing energy levels have been esti- 
mated by Drellishak, et al. (ref. 8) and have been published (ref. 2) up to a principal 
quantum number of 32. 


A survey of methods for treating the excitation partition function and determining 
the principal quantum number cutoff is presented in reference 16. Based on the theory of 
Lewis and Margenau (which was also used in ref. 8) the principal quantum number cut- 
off n is given by 


n « 30(q + 1) 


1/2 


10 


-I V4 


n c 


U0 
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(A2) 


where q = 0 for an atom and q = 1 for an ion. The equation is made dimensionless 
for easy reference so that at T = 10^ °K and n e = 10^2 m“3, n ~ 30(q + 1)1/^. 

Equation (A2) is valid as written only for a mixture of neutrals and equal numbers 
of electrons and single ions. Shown in figure 30 are values of n for the temperature- 
pressure region considered and degrees of seeding of K c = 0.01 (X = Cesium) and 
K c = 0. 
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For the alkali and nitrogen atoms the evaluation of the electronic excitation parti- 
tion function as written in equation (8b) requires, particularly at large values of n, a 
large amount of input data (for example, at n = 20, for cesium N(n) = 199) and a corre- 
spondingly long computational time. Because of the functional dependence of Qg(i) 
and n, the determination of n is accomplished with an iteration procedure (see the 
section "Determination of Species Concentrations”) and it is extremely important that 
the time required to compute Qg(i) be reduced to a minimum. With this in mind it is 
advantageous to write 

N(n) _ 

Q e (i) = £ gj, e T = Q e (0 1 + Q e (i) 2 (A3) 

v=l 

where 

N ( n o) _ 2a 

Q e (i) l s Yj gpQ T 

v=l 

and is a function of temperature only, and where 

N(n) n 

Qe(i) 2 = ^ T = ^ a(2m2)e 

te=N(n 0 )+l m=n 0 +l 

and is a function of both temperature and n. The quantity represents either 0 2 (x) 
or 03 (N) and a is the number of parent configurations (a = 1 for the alkali-metal 
atoms and a = 9 for the nitrogen atoms). The quantity Qe(i)j is a summation over 
all energy states with a principal quantum number m such that m g n 0 and Q e (i) 2 
is the contribution to Qg(i) resulting from a summation over the states where 
n D < m s n. The principal quantum number n 0 is chosen so that all states with m > n 0 
will be approximately hydrogen-like and also so that for a particular m the energy 
states for each azimuthal quantum number will differ by very small amounts of energy. 
The total contribution to Q^i^ for each principal quantum number greater than n Q is 
then approximately 

* 9 H 

_ * m2 

2m2e T 

where 2m2 is the sum of degeneracies for all azimuthal quantum numbers. 



24 



APPENDIX A 


A simple expression may be obtained for Qg(i) 2 if, in addition, n Q is chosen 
so that 



« 1 


and A0 m is defined so that 



where n s is some reference term n Q < n s < n. The quantity Q^i^ may now be 
written as 

_ 0s n _ 

Qe(i) 2 = 2ae T ^ m2 e T 
m=n 0 +l 


Since 


0s 


e 


0 H 






1 


Q e (i) 2 may be written as 

Oe(i) 2 “ 2ae T 

m=n 0 +l 


If the reference term n s is chosen (thereby defining A0 m ) so that 
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-and, if the definition of A0 m is recalled, 

a a 


n 


n 


I 


m 


2 _ 


■ I 


m 


m=n 0 +l 


m=n 0 +l 


2 1 _ _1 
T (n g 2 m2 


= 0 


Then n g ^ becomes 


I 




_ m=n 0 +l _ 2 (n 2 + nn Q + n 2 ) - 3(n + n 0 ) + 1 

6 


n s 2 = 

s n - n, 


and Qe(i )2 becomes 


T 


Q e (i) 2 - 2ae (n - n 0 )n g ' 


(A4) 


Equation (A4) provides a simple and accurate manner of estimating the contribution to 
Q e (i) caused by higher quantum numbers. For both the nitrogen and alkali-metal atoms, 
n D was taken as 15 although for the higher temperatures considered, an n Q less than 15 
could have been used with sufficient accuracy for the alkali metals. Shown in table IV 
are the values of g v and 8 V used to evaluate Q e (i)j- In many cases approximately 
equal energy states have been combined to reduce the amount of input data. The data are 
based on all energy levels with a principal quantum number less than or equal to 15. For 
the ions it is sufficiently accurate to take Q^ijg = 0. 

The quantities e e (i) and cr e (i) can be obtained in a manner similar to that used 
to obtain equations (A3) and (A4) and they are 



and 


o-e(i) = 



- e e 2 (i) 


(A6) 
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The excitation of species N2 becomes significant only in temperature-pressure regions 
where n^ « nt and therefore over the whole temperature-pressure region, it was 
sufficiently accurate to assume that 


Qe( N 2) = 1 

% (N 2 ) = <J e (N 2 ) = 0 

By chain rule differentiation it is easily established that for a general function 
f(T,n) of temperature and principal quantum number cutoff n, 



where n = Constant (T/n e )*/^ and is treated as a continuous, rather than discrete, 
variable. For the three functions Q e (i), eeU)> and a e (i) which are of the form f(T,n), 
the second term in equation (AT) is generally much smaller than the first. Therefore, at 
a given temperature and n 

V /p or p 


'9f(T,n)\ 

. 3T 
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SIMPLIFICATION OF THE NUMBER DENSITY EQUATIONS 


The number density equations written in terms of a r are given by equations (14a) 
to (14c). The equations in this form contain none of the usual assumptions of indepen- 
dence of reaction. (See, for example, ref. 3.) By combining equations (14a) and (14c), 
it is possible to obtain 


ao = 


1 + 


/Ki(N)/l - aj) 


K 3 (N)\ 2a 


1 / Ur 


1-Q! 


i 1 - a i) 


K C 1 
— + 1 

“3 


(Bl) 


and in a similar manner equations (14b) and (14c) together produce 

a 1 

3 ■ | KaPO fi - 
K 3 (nA a 2 I 

Equations (Bl) and (B2), respectively, provide the useful inequalities 



and 


“3 




K 3 (N) 

k 2 (x) 


(B2) 


(B3) 


(B4) 


In addition, an inspection of equations (14a) to (14c) reveals that for K r (Y) » P, the 
corresponding reaction r is essentially complete (a r ~ l). Therefore from equa- 
tions (B3), (B4), and figure 11, it follows that in the regions of reaction (0 < a r < l), 


“st 1 ' a 2) « 1 

o^ 1 - a x ) « 1 


which is true only, if for 0 < aj < 1 
then Q!j ~ 1 and ~ 1. 


and 0 < 0*2 < 1, then ag 


0 or, if for 0 < ar 3 < 1 
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Although it may be expected that reaction 2 will become significant at lower tem- 
peratures than reaction 1, this is not necessarily so. By using conditions (B5), equa- 
tions (14a) and (14b) may be combined to provide 


1 - a\ _ 2K 2 (X) 1 - 02 
t t 2 K c Kj(N) a 2 2 


(B5) 


Q!i 


This relation in connection with figure 11 illustrates that, particularly at larger values 
of K c , reaction 1 may become significant in regions where reaction 2 is not essentially 
complete. (See fig. 1.) 

Based on this discussion, equations (15) to (17) were obtained from the following 
equations: 


aj2 

_ K X (N) 

£ + a l 

+ 2K C (1 + a 2 ) 

] (B6) 

1 " a l 

4P 

a 2 2 

k 2 (x). 

1 + 

+ 2K c (1 hh a 2 ) 

(B7) 

l-oi2' 

" 2K c P 1 

a 3 (a 3 

1 - 

+ K c)_ 
a 3 

K 3 (N) 

P 

f + “3 + 2K c) 

(B8) 


The quantities <p j to 0 5 are given by 

, “it 1 ~ a i l 

2 - oi-i 


_ Oi 2 (l ~ 0! 2 ) 

^2 = ^^ 


a 3( a 3 + KcX 1 - + a 3 + 2K c) 

(-K c a 3 2 + 4Kca 3 + 2a 3 + K c + 2K C 2 ) 

_ [l + “1 + 2K c(! +a 2i| ( 2 - a l)( 2 - «2) 


04 = 


2 [2 - Q-2 +2K C (2 - o^)] 


a a ia 2 (l ~ «l )(l " «2) 

5 2[2 - a 2 +2K C (2 - a{( 


29 



APPENDIX C 


COMPUTER PROGRAM OUTLINE 

As mentioned in the text, many high-enthalpy gas flow problems require, for a 
solution, a high-speed computer program using the thermodynamic properties of the gas 
in connection with a step-by-step integration procedure. This outline is presented in 
the hope that it might be useful either as a subroutine to supply all necessary thermo- 
dynamic variables at each step in the integration or for preparing tables of thermody- 
namic properties. 

(1) Read in fundamental and gas dependent constants, including 7r, k, h, m n , 

me> Pq> -rot’ @y> ®r(^> ^H’ ^w* 

(2) Read in excitation constants g v and 6 V for v = 1 to v = N(n 0 ) for each 
species. (See table IV for values of g v and 9 V .) 

(3) Enter subroutine with values of K c , T, and P/= £-\. 

N ( n o) b _0v ' ° 

(4) Compute ^ gj,^j e T ( b = 2 ) for each species. 

(5) Start the iteration on n by choosing an initial n = n^ = 20) and computing 
Q e (i) from equation (A3), and K r (Y). (See section "Equilibrium Constants.") 

(6) Start iteration on A by choosing A = Aj (Aj = 2Kc) and computing ai 
(eq. (15)) and a 2 (eq. (16)). 

(7) Compute a|= 2K c | H A = Af, go to step 8; if A + Aj, set Ai = A 

and go to step 6. ' ^ 

(8) Compute a% (eq. (17)), Z (eq. (19)), and n e (eq. (18)). 

(9) Compute n (eq. (A2)). If n = n 1? go to step 10; if n * ni set nj = n and 
go to step 5. 

(10) With the correct values of a r and n, compute e e (i) (from eq. (A5)), a e (i) 
(from eq. (A6)), e v , cr v (see section "Partition Functions"), e ro t, and a ro j-. 

(11) Compute R (eq. (6a)), j3 1>2j3 (eqs. (12a), (12b), and (12c)), ju 12>3 

(eqs. (13a), (13b), and (13c)), (eqs. (20a), (20c), and (20e)), (eqs. ( 20b )> 

(20d), and (20f)), /3 Q , and p 0 (see definition). 
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(12) Compute any of the desired properties from among the following: p (eq. (5)); 
nj (eq. (18)); H, E, and S (see section "Internal Energy, Enthalpy, and Entropy"); 

Cp, Z T , Zp, y, a, and 0 (see section "Specific Heats and Speed of Sound"). 
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TABLE I.- CHARACTERISTIC TEMPERATURES FOR DISSOCIATION AND IONIZATION 

[fl rot = 5.784° K; 0 V = 3395° k] 


Gas, Y 

r ~ 

K w 

Characteristic temperature of — 

Dissociation, 

°K 

First 

ionization, 

°K 

Second 

ionization, 

°K 

Nitrogen, N 


113 220 

168 730 

343 560 

Lithium, Li 

0.4954 

11 950 

62 550 

877 550 

Sodium, Na 

1.641 

8 470 

59 630 

548 790 

Potassium, K 

2.791 

5 965 

50 350 

369 150 

Rubidium, Rb 

6.102 

5 690 

48 460 

319 130 

Cesium, Cs 

9.488 

5 220 

45 180 

291 280 
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TABLE n.- COMPARISON OF PRESENT RESULTS WITH REFERENCE RESULTS 



T,°K 
2 000 

Present 

paper 

1.000 

Z 

Reference 2 
1.000 

Reference 3 
1.000 

Present 

paper 

P = 

0.2303 X 10 7 

H, J/kg 

Reference 2 Reference 3 
0.01 atmosphere 
0.2314 X 10 7 | 0.2302 X 10 7 

Present 

paper 

1.273 x 10 3 

c p , J/kg-°K 
Reference 2 

1.281 X 103 


4 000 

1.009 

1.009 

1.009 

.5220 

.5257 

.5222 

2.461 

2.507 


6 000 

1.774 

1.765 

1.773 

3.494 

3.468 

3.489 

20.08 

22.31 


8 000 

2.045 

2.045 

2.045 

4.940 

4.939 

4.932 

6.158 

6.194 

* 

10 000 

2.481 

2.484 

2.482 

7.885 

7.864 

7.838 

27.43 

27.08 


12 000 

3.530 

3.546 

3.554 

14.52 

14.52 

14.49 

26.33 

26.66 


14 000 

3.904 

3.927 

3.939 

17.30 

17.33 

17.34 

6.509 

6.571 

t 

16 000 

3.958 

3.982 

3.992 

18.23 

18.29 

18.28 

3.711 

4.123 


2 000 

l.ooo 1 

1.000 

1.000 

P = 

0.2303 X 10 7 

1.0 atmospher 
0.2314 x 10 7 

e 

0.2302 x 10 7 

1.273 x 10 3 

1.281 x 10 3 


4 000 

1.001 

1.001 

1.001 

.4938 

.4981 

.4937 

1.432 

1.457 


6 000 

1.121 

1.118 

1.121 

1.183 

1.180 

1.182 

8.364 

8.162 


8 000 

1.839 ; 

1.830 

1.839 

4.104 

4.076 

4.092 

11.30 

11.63 


10 000 

2.040 

2.040 

2.040 

5.367 

5.369 

5.343 

5.494 

5.483 


12 000 

2.243 

2.246 

2.244 

7.005 

7.027 

6.953 

12.18 

12.18 


14 000 

2.727 

2.735 

2.735 

10.40 

10.47 

10.28 

21.59 

21.72 


16 000 

3.353 

3.352 

3.385 

14.74 

14.82 

14.67 

19.20 

19.00 


18 000 

3.723 

3.711 

3.778 

17.59 

17.62 

17.63 

9.866 

9.626 


20 000 

3.861 

3.849 

3.923 

19.04 

19.10 

19.14 

5.512 

5.847 


Reference 3 


1.273 X 10 3 
2.459 
19.90 
5.741 
26.92 
26.83 
6.583 
3.670 


1.273 x 10 3 
1.433 
8.346 
10.69 
4.524 
11.62 
21.28 
19.71 
10.21 
5.518 
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TABLE HI.- COMPARISON OF DATA BASED ON BOTH A FIXED AND A VARIABLE PRINCIPAL 

QUANTUM NUMBER CUTOFF 


(a) Unseeded 


T, °K 

ne, m-3 


Z 

H, J/kg 

Cp, J/kg-°K 


L n = 5 

Variable n 

n = 5 

Variable n 

n = 5 

Variable n 

n = 5 

Variable n 





P = 0.01 atmosphere 




o nrvn 



1 000 

1 000 

0 2303 X 10 7 

0.2303 X 10 7 

1.273 x 10 3 

1.273 X 10 3 

a nnn 



1 009 

1 009 

5220 

.5220 

2.461 

2.461 

6 000 

5.276 X1018 

5.276 X 10 18 

1.774 

1.774 

3.494 

3.494 

20.08 

20.08 

8 000 

.2099 X 1021 

.2099 X 1021 

2.045 

2.045 

4.936 

4.940 

6.061 

6.158 

10 000 

1.427 

1.423 

2.483 

2.481 

7.841 

7.885 

26.94 

27.43 

12 000 

2.674 

2.651 

3.554 

3.530 

14.49 

14.52 

26.82 

26.33 

14 000 

2.581 

2.557 

3.939 

3.904 

17.34 

17.30 

6.585 

6.509 

16 000 

2.287 

2.269 

3.989 

3.958 

18.27 

18.23 

3.686 

3.711 





P = 1.0 atmosphere 




2 000 



1.000 

1 000 

0.2303 X 10 7 

0.2303 X 10 7 

1.273 X 10 3 

1.273 X 10 3 

4 000 



1.001 

1 001 

.4938 

.4938 

1.432 

1.432 

6 000 



1.121 

1.121 

1.183 

1.183 

8.364 

8.364 

8 000 

0.1953 X 10 22 

0.1953 X 10 22 

1.839 

1.839 

4.103 

4.104 

11.28 

11.30 

10 000 

1.772 

1.771 

2.040 

2.040 

5.357 

5.367 

5.333 

5.494 

12 000 

6.671 

6.652 

2.244 

2.243 

6.955 

7.005 

11.71 

12.18 

14 000 

14.09 

13.97 

2.735 

2.727 

10.29 

10.40 

21.31 

21.59 

16 000 

18.76 

18.51 

3.384 

3.353 

14.68 

14.74 

19.69 

19.20 

18 000 

19.18 

18.87 

3.776 

3.723 

17.62 

17.59 

10.20 

9.866 

20 000 

17.97 

17.69 

3.918 

3.861 

19.12 

19.04 

5.585 

5.512 











(b) 5-pereent seeded 


T, °K 

Re* 

m“ 3 


Z 

H, J/kg 

Cp, j/kg-°K 

n = 5 

| Variable n 

n = 5 

Variable n 

n = 5 

Variable n 

1 

j n = 5 

j Variable n 





P = 0.01 atmosphere 




2 000 

0.1054 x 10 20 

0.1054 X 1020 

1,100 

1.100 

0.1666 x 10 7 

0.1666 X 10 7 

0.9342 X 10 3 

0.9359 X 10 3 

4 000 

13.05 

12.55 

1.195 

1.191 

.4706 

.4705 

2.259 

2.232 

6 000 

6.144 

5.982 

1.989 

1.986 

2.556 

2.555 

13.01 

13.02 

8 000 

4.966 

4.896 

2.218 

2.216 

3.419 

3.420 

3.568 

3.605 

10 000 

15.37 

15.32 

2.656 

2.654 

5.408 

5.436 

18.98 

19.29 

12 000 

26.93 

26.71 

3.752 

3.728 

10.03 

10.05 

18.49 

18.15 

14 000 

25.83 

25.60 

4.139 

4.104 

11.99 

11.97 

4.575 

4.524 

16 000 

22.88 

22.70 

4.189 

4.157 

12.64 

12.61 

2.608 

2.625 





P = 1.0 atmosphere 




2 000 

0.0106 X 1022 

0.0106 X 1022 

1.100 

1.100 

0.1663 X 10 7 

0.1663 X 10 7 

0.9161 x 10 3 

0.9172 X 10 3 

4 000 

3.129 

3.020 

1,120 

1.119 

.3774 

.3801 

1.373 

1.407 

6 000 

8.015 

7.252 

1.317 

1.307 

.9606 

.9576 

6.315 

6.273 

8 000 

4.442 

4.144 

2.046 

2.039 

2.968 

2.964 

7.166 

7.190 

10 000 

4.038 

3.889 

2.214 

2.209 

3.721 

3.722 

3.277 

! 3.341 

12 000 

7.889 

7.825 

2.410 

2.407 

4.797 

4.825 

8.240 

8.524 

14 000 

14.67 

14.55 

2.916 

2.907 

7.138 

7.205 

14.95 

15.14 

16 000 

18.96 

18.72 

3.580 

3.547 

10.19 

10.23 

13.64 

13.31 

18 000 

19.23 

18.94 

3.975 

3.922 

12.23 

12.20 

7.069 

6.842 

20 000 

17.98 

17.72 

4.118 

4.060 

13.27 

13.21 

3.926 

3.878 
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TABLE IV.- CHARACTERISTIC TEMPERATURES FOR ELECTRONIC EXCITATION 


(a) Atoms 


Index 

Nitr 

ogen, N 

Ces 

ium, Cs 

Potassium, K 

Sodium, Na 

V 

g* 

e v 

Si/ 

K 


9 v 

gy 

*v 

1 

4 

0 

2 

0 

2 

0 

2 

0 

2 

10 

27 662 

6 

16 614 

6 

18 739 

6 

24 414 

3 

6 

41 495 

10 

20 946 

12 

30 862 

2 

37 035 

4 

12 

119 910 

2 

26 670 

6 

35 559 

10 

41 974 

5 

6 

124 010 

6 

31 491 

12 

39 433 

6 

43 554 

6 

12 

126 810 

10 

32 538 

14 

40 471 

2 

47 770 

7 

22 

136 290 

16 

35 182 

6 

41 732 

30 

49 877 

8 

12 

137 420 

16 

37 335 

44 

43 885 

2 

52 333 

9 

14 

139 290 

34 

38 830 

6 

44 708 

48 

53 356 

10 

6 

140 710 

16 

39 941 

66 

45 916 

66 

55 233 

11 

10 

143 390 

62 

40 837 

18 

46 734 

8 

55 607 

12 

22 

149 210 

90 

41 909 

98 

47 237 

98 

56 434 

13 

86 

150 910 

120 

42 689 

270 

48 200 

288 

57 464 

14 

22 

153 650 

326 

43 445 

192 

48 790 

192 

58 063 

15 

32 

154 990 

242 

43 793 

242 

49 046 

242 

58 325 

16 

104 

158 420 

288 

44 036 

288 

49 263 

288 

58 536 

17 

54 

160 380 

338 

44 219 

338 

49 430 

338 

58 700 

18 

40 

162 100 

392 

44 358 

392 

49 561 

392 

58 830 

19 

416 

162 670 

450 

44 470 

450 

49 666 

450 

58 935 

20 

608 

164 590 







21 

966 

165 700 







22 

2552 

166 710 







23 

3366 

167 200 







24 

6570 

167 560 











(b) Ions 





Index 

Nitrogen, N + 

Cesium, Cs + 

Potassium, K* 

Sodium, Na + 

V 



Sy 


gy 

e v 

gy 

e v 

1 

1 

0 

1 

0 

1 

0 

1 

0 

2 

3 

70.65 

10 

152 560 

9 

234 310 

9 

381 770 

3 

5 

188.9 

14 

155 160 

12 

237 200 

3 

386 700 

4 

5 

22 036 

3 

159 630 

3 

239 510 

3 

421 890 

5 

1 

47 030 

14 

161 770 

16 

244 650 

15 

428 140 

6 

5 

67 865 

7 

163 620 

12 

247 190 

9 

430 480 

7 

15 

132 720 

10 

172 710 

3 

263 600 

4 

432 330 

8 

9 

157 140 

4 

176 600 

12 

268 360 

1 

444 400 

9 

5 

207 460 

3 

177 890 

8 

270 380 

25 

459 120 

10 

12 

214 460 

3 

182 040 

12 

272 920 

23 

477 320 

11 

3 

223 200 

12 

185 150 

1 

280 240 

24 

479 020 

12 

3 

236 840 

15 

187 410 

8 

306 080 

15 

508 620 

13 

21 

240 210 

1 

191 580 

13 

309 940 



14 

9 

245 510 

3 

203 670 

8 

311 740 



15 

5 

250 660 

9 

206 480 

16 

312 710 



16 

1 

256 500 

12 

214 390 

12 

316 050 



17 

41 

269 020 

27 

220 490 

8 

320 920 
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Figure 9.- The electronic excitation partition function derivative o e (i) for the nitrogen atoms. 
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